Reynolds number effect on the velocity increment skewness in isotropic turbulence 
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Second and third order longitudinal structure functions and wavenumber spectra of isotropic 
turbulence are computed using the EDQNM model and compared to results of the multifractal 
formalism. At the highest Reynolds number available in windtunnel experiments, R\ — 2500, both 
the multifractal model and EDQNM give power-law corrections to the inertial range scaling of 
the velocity increment skewness. For EDQNM, this correction is a finite Reynolds number effect, 
whereas for the multifractal formalism it is an intermittency correction that persists at any high 
Reynolds number. Furthermore, the two approaches yield realistic behavior of second and third order 
statistics of the velocity fluctuations in the dissipative and near-dissipative ranges. Similarities and 
differences are highlighted, in particular the Reynolds number dependence. 



I. INTRODUCTION 

The nonlinearity in the Navier-Stokes equations gives 
rise to an interaction between different length-scales in a 
turbulent flow. These interactions are the basic mecha- 
nism behind the celebrated Kolmogorov-Richardson en- 
ergy cascade [l|, This phenomenological picture of 
energy cascading from scale to scale towards the scales 
in which dissipation becomes appreciable is the corner- 
stone of a large number of turbulence models (e.g. ref- 
erence 0). If locality in scale-space is assumed, energy- 
conservation and local-isotropy will lead to a wavenum- 
ber dependence of the energy spectrum of the form 



,2/3^„-5/3 



(1) 



with e the energy flux, which, using the assumption of 
statistical stationarity, equals the energy dissipation. A 
physical space equivalent of this scaling law is the scale 
dependence of the second-order longitudinal structure 
function, 



DLLir) 
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The definitions of DLL{r) and E{k) will be given below. 

The possibility of corrections to the inertial range scal- 
ing of structure functions, due to the intermittent char- 
acter of the energy dissipation jj, was taken into ac- 
count in a more general theory advanced by Kolmogorov 
and Oboukhov p, Q . Experiments aiming at the mea- 
surement of the intermittency corrections (e.g. reference 
0, Si) indeed showed small corrections to the scaling 
which could be due to intermittency, in particular for 
higher-order structure functions. Subsequently a large 
number of phenomenological models was proposed to de- 
scribe the intermittent character of turbulence. Refer- 
ence 0] gives an overview of work on intermittency upto 
1995. One of the more successful models, in the sense 
of reproducing the different features of isotropic turbu- 
lence, is the multifractal model Q . This phenomenologi- 
cal description compares well to measurements and gives 



non-zero intermittency corrections to the inertial range 
scaling of the energy spectrum and of higher order quan- 
tities. 

A valuable theoretical tool to study the statistical 
properties of homogeneous turbulence is two-point clo- 
sure theory. The first theoretical approach of this kind, 
derived from the Navier-Stokes equations, is the Direct 
Interaction Approximation (DIA) [l3|- Subsequent im- 
provements [ll| of this theory allowed to show that the 
^-5/3 dependence of the energy spectrum can be related 
directly to the Navier-Stokes equations. Simplifications 
led to different related closures such as the test-field 
model [l^ . the Lagrangian renormalized approximation 
p!3 | and the Eddy-Damped Quasi-Normal Approxima- 
tion (EDQNM) [14,]. EDQNM is of the closures named 
here the simplest. It is obtained by assuming in the DIA 
formulation that the two-time correlations decay expo- 
nentially with a typical time-scale modeled phenomeno- 
logically. We note that this time-scale can also be de- 
termined self-consistently within the EDQNM approach 

These closures, although directly related to the Navier- 
Stokes equations, do not yield any intermittency correc- 
tions to the scaling ([T]). However, the predicted results for 
scaling exponents of the energy spectrum compare rather 
well to experimentally observed values [l^. Indeed, at 
low Reynolds numbers, corrections to the scaling expo- 
nents due to the finite Reynolds number are usually larger 
than the expected intermittency corrections and these 
finite Reynolds number effects vanish very slowly [l6l - 
fisj . The fact that two-point closure and the multifractal 
formalism can treat both low and very high Reynolds 
numbers using limited computational effort, makes these 
approaches very attractive to study Reynolds-number ef- 
fects. 

The present work will compare the predictions of clo- 
sure for second and third-order quantities with results 
of the multifractal description. This will allow to show 
to what extent intermittency corrections can be distin- 
guished from Reynolds number effects at low, moder- 
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ate and high Reynolds numbers. We choose to compare 
quantities in physical space, since most experimental and 
theoretical efforts aiming at the understanding and de- 
scription of intermittency focus on these quantities (we 
note however that in principle intermittency corrections, 
if any, should also be observed in wavenumber spectra). 
Therefore we need to convert the Fourier-space quanti- 
ties into physical space quantities. It is described in the 
next section how this is done. The relations to convert 
physical space quantities into their Fourier-space coun- 
terparts is also given. In section Hill we will present the 
EDQNM model and we will give an outline of the multi- 
fractal description. In section HVl we present the results 
of the EDQNM model for these quantities and compare 
with the multifractal results in both Fourier and physical 
space. Section IVl concludes this article. 



II. EXACT RELATIONS BETWEEN SECOND 

AND THIRD ORDER QUANTITIES IN 
FOURIER SPACE AND IN PHYSICAL SPACE 



In this section we will give the relation between the 
energy spectrum and the second order structure func- 
tion Z?ii(r), and between the nonlinear transfer and 
the third-order longitudinal velocity structure function, 
DLLL{f). Even though the relations given here are not 
new (e.g. d^-HH), the details of the derivation are dis- 
persed or not well documented in literature and we think 
that it is therefore worth to write down in detail this 
derivation, which can be found in the appendix. 



A. Derivation of the Lin-equation 

The starting point is the Navier-Stokes equations for 
incompressible flow. 



with 



dui{x) 
dt 



■Uj{x) 



dui(x) 1 dp{x) d'^Ui{x) 
dxj p dxi 



(3) 



dujjx) 
dxi 



= 0, (4) 



with p the density and p the pressure. Time arguments 
are omitted for brevity. The three-dimensional Fourier 
transfer is defined as 



u^{k) ^ j^-^ j u,{x)e rffc. 



(5) 



In Fourier space the Navier-Stokes equations can be writ- 
ten 

dui{k) 2 



dt 



+ vk Ui{k) 



'-2^i3m{k) If Uj{p)umiq)S{k - p - q)dpdq (6) 



Pijm{k) — kjPi„i(k) + kj^Pij^k), 



Pij {k^ — ^i' 



ki kj 



(7) 
(8) 



To derive ([B]), the incompressibility condition was used 
to eliminate the pressure term. In isotropic non-helical 
turbulence, the energy spectrum is related to Uiik) by 



P,,{k) 



E{k) = u,{k)u*(k) 



and since Ui{x) is real, this gives 

E{k) = 27rfc2u,(fc)u,(-fc). 



(9) 



(10) 



In order to derive the equation for E(k), we multiply 
([6]) by Uil—k). Then we write a similar equation for 
Ui{~k) and multiply by Ui{k) Summing both equations 
and averaging yields. 



1 + 2^^' 
at 



E{k) = l^eP,,m{k)[T,,ra{k)~T:^^{k)\ 

= T{k) (11) 



with 



(k) = JJ Ui{k)uj{p)u„i{q)S{k + p + q)dpdq 



(12) 



Tum{k)= // Ui{-k)uj{-p)u„,[-q)5{k + p + q)dpdq. 



(13) 



By isotropy it can be shown that T*^^{k) = —Tij.^[k). 
The RHS of the Lin-equation ([TT|) is the nonlinear trans- 
fer T{k), which we will relate to the third-order longitudi- 
nal structure function. But first we will give the relation 
between the energy spectrum and the second order lon- 
gitudinal structure function. 



B. Relation between E{k) and Dll{t). 

The second order longitudinal structure function is de- 
fined as 



DLL{r) = 5ul 



(14) 



with 



duL ^ UL - u'l^ = —Ui{x) -Ui{x + r). (15) 



its relation to the energy spectrum is 



DLLir) 



E{k)f{kr)dk 



(16) 
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with f{x) given by 



1 sin(a;) — (x) cos(x) 
3 (x)^ ' 



(17) 



Using this in (ITBl) and ([20)) . we find for very smaU r, 
2 



DLLir) = —r' I k'E{k)dk 



The derivation of this expression is given in the appendix. 
A convenient expression to compute the energy spectrum 
from the second order structure function is 



er 
15^ 



DLLL{r) - j T{k)Ak~ ^r'' j k^T{k)dk 



1 1 kr [sin(fcr) — fcr cos(fcr)] dr 



2m2 



35 



k'^T(k)dk 



(18) in which we used that 



and we refer to Mathieu and Scott [211 for the derivation. 



C. Relation between T{k) and Dlll{j). 

The third order longitudinal structure function in ho- 
mogeneous turbulence can be expressed as 



DLLL{r) = 5u\ = 3 [ulu'j^ - u'j^u\ 
which is related to the transfer spectrum by 



(19) 



DLLL{r) = r / T{k)g{kr)Ak, 
Jo 



with 



g{x) = 12 



3 (sin X — X cos x) — x'^ sin x 



(21) 



with details given in the appendix. The equivalent ex- 
pression to compute the transfer spectrum from DLLL{r) 
is 1211, 



T(fc) 



6n 



sin(A:r) d 
r dr 



ld_ 

r dr 



[r'DLLLir)) 



dr. (22) 



D. Small scale behavior of DLL^r) and DLLL(r) 

Before continuing, let us have a look at the behavior 
of the functions ^} and (PT|). 



f{x) = 4 



1 sin X — X cos X 
3 ^3 



(23) 



, , „ 3 (sin a; — x cos x) — sin a: 
g{x) = 12^ ^ . (24) 



Taylor expansions of the sine and cosine terms show that 
for a; J, 0, 



f{x) = -x^ + 0{x^) 
g{x) = ^~^x^ + 0{x''). 



(25) 
(26) 



2v j k^E{k)Ak = e 
j T{k)dk = 0, 



(27) 

(28) 

(29) 
(30) 



with e the energy dissipation. So we find that the struc- 
ture functions of order 2 and 3 scale as and respec- 
tively for very small r, which is expected since at small 
enough scales the flow can be considered as smooth. 
The velocity-increment skewness is defined as 



S{r) 



DLLLjr) 



(20) Since at very small scales 



du 

duL ~ r— , 
ox 



one finds that 



lim S(r) 



{du/dxY 



3/2 ■ 



(31) 



(32) 



(33) 



{du/dxYj 

Using expressions (|27| and ([28)) . we find [l| 



lim S{r) 



153/2 eT{k)dk 



r-!-0 



P,u)2]3/2 35(2)1/2 [J- k^E{k)dkY' 



P4) 



In the case of high-Reynolds number, if the non- 
stationarity can be neglected at high fc, or if the tur- 
bulence is kept stationary by a forcing term acting only 
at small fc, we have 



j k'^T{k)Ak « j 2vk'^E{k)Ak, 



(35) 



so that 191 



153/221/2^ j^k^E{k)dk 
^To^W« 35 [J-k^E{k)dk]m - (36) 

The velocity-derivative skewness is then completely de- 
termined by moments of the energy spectrum. Using 
(j29l) , and assuming an inertial range spectrum extending 
upto fc/ of the form 



E{k) ~ e2/3r 



(37) 
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we obtain for a > —3 



l/(a+3) 



(38) 



Substituting this in p6p we obtain, 



lim S{r) 

r—^O 



-|(5+3a) 



3q + 5 



(39) 



in which we used that the Taylor-scale Reynolds number, 
to be defined later, is proportional to v^^^^. It follows 
from this expression that the velocity derivative skewness 
is independent of the Reynolds number if a = —5/3. If 
corrections to the Kolmogorov 1941 (K41) scaling are 
present, as is the case in the multi- fractal formalism, the 
skewness becomes a function of the Reynolds number. 
This dependence is by (|39l) directly related to the inter- 
mittency correction to the K41 scaling. 



E. Large scale behavior of DLL{r) and Dlll{t) 



At large r we find 
DLL{r) 



E{k)dk 



(40) 



which is expected from (lAip since at large separation dis- 
tances the correlation between the velocity at two points 
is supposed to vanish. DLLL{r) tends for the same rea- 
son to zero for large separation distances r. At large 
r the velocity increment skewness should therefore go 
smoothly to zero, since the velocity correlation should de- 
cay smoothly at large r. The exact way in which DLLL{r) 
tends to zero depends on the behavior of the energy spec- 
trum at the very low wavenumbers. 



III. THE EDQNM MODEL AND THE 
MULTIFRACTAL DESCRIPTION 

A. The EDQNM model 

The EDQNM model is a closure of the Lin-equation 
in which the nonlinear transfer T[k) is expressed as a 
function of the energy spectrum. The transfer T{k) is 
given by 

T{k) = ff Qkp, {xy + z^) [k^p E{p)E{q) 



A 



^p'E{q)Eik)] 



dpdq 
pq 



(41) 



In equation (|41l) , A is a band in p, g-space so that the 
three wave- vectors fe,p, q form a triangle, x, y, z are the 
cosines of the angles opposite to k, p, q in this triangle. 
This particular structure is common to all closures de- 
rived from the Direct Interaction Approximation [lo| . 



DIA is a self-consistent two-point two-time theory with- 
out adjustable parameters. Simplifications are needed 
to obtain a single-time (or Markovian) description intro- 
ducting assumptions and adjustable parameters. In the 
case of EDQNM the simplifying assumption is that all 
time-correlations decay exponentially, with a time-scale 
Qkpq modeled phenomenologically by 



e 



kpq 



1 - exp(-(77fc + r]p + r]q) x t) 

Ilk + lip + Vq 

in which rj is the eddy damping, expressed as 



(42) 




'-E{s)ds + iyk'^ 



(43) 



related to the timescale associated with an eddy at 
wavenumber k, parameterised by the EDQNM parame- 
ter, A, which is chosen equal to 0.49 [l5|. The exponential 
time-dependence in (j42[) appears by the assumption that 
the initial conditions have vanishing triple correlation as 
would be the case for a Gaussian field. Its influence van- 
ishes at long times. For an extensive discussion of the 
EDQNM model see [22l, '23'], but we want to stress that 
one of the key features of EDQNM is that it is applicable 
at all Reynolds numbers (it is not an asymptotic theory) 
and at all scales of a turbulent flow. In other words, its 
results go beyond mere scaling and can give insights on 
the Reynolds number dependency of different quantities 
related to turbulence. 

We performed simulations of the EDQNM model in 
the unforced case by integrating numerically Eq. (j4ip . 
starting from an initial spectrum. 



Ek{0) = BkUxp [-{k/kLf], 



(44) 



with B chosen to normalize the energy to unity and k^ — 
Sko, ko being the smallest wavenumber. The resolution is 
chosen 12 gridpoints per decade, logarithmically spaced. 
In the decaying simulations results are evaluated in the 
self-similar stage of decay, in which e/ekin, with ekin the 
kinetic energy, is proportional to t^^. Forced simulations 
are evaluated when a steady state is reached. The forcing 
in these cases corresponds to a region k < k^, in which 
the energy spectrum is kept constant. 



B. The multifractal description 

In this section some key concepts of the multi-fractal 
description will be presented. A more detailed presenta- 
tion and references can be found in appendix [Ul 

In the multifractal description the velocity-increments 
Sul{x, r) = ul{x + r) — ul{x) at scales r in the inertial 
range are modeled by the product of two independent 
random variables 



(45) 
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In this expression ^ is a zero average Gaussian random 
variable of variance , where is twice the mean-square 
of the velocity fluctuations. The quantity (3^ introduces 
the scale dependence in the statistics of the velocity in- 
crements. It is defined as 



Pr = 



(46) 



with Lq the integral lengthscale. The particularity of the 
approach lies in the fact that the exponent /i is a fluc- 
tuating quantity. If a constant value /i = 1/3 is chosen, 
K41 behavior is recovered. In the multifractal framework 
h is determined by the probability density function 



Vrih) <x — 

>0 



l-V(h) 



(47) 



If the unknown function I?(/i) is given (for example by 
comparison with experimental results), a complete de- 
scription of the inertial range statistics of the velocity 
increments can be obtained. An extension to take into 
account the dissipative effects was proposed by Paladin 
and Vulpiani 2j] and Nelkin 25| . Further details on the 
multi-fractal description, including the expressions for /3r 
and Prih), are given in the appendix. 



IV. RESULTS FOR SECOND AND THIRD 
ORDER QUANTITIES 

A. Results in Fourier space 

In the following we will consider three different values 
of the Reynolds number 



15- 



(48) 



These values are Rx = 380, 2500, 25000. The low- 
est corresponds to a typical Reynolds number for lab- 
oratory experiments in jets or wind-tunnels, the second 
one to the highest Reynolds number obtained in wind- 
tunnel turbulence, i.e. in the Modane windtunnel j26| . 
and R\ = 25000 corresponds to the Reynolds number of 
large scale atmospheric flows and no controlled experi- 
mental results of isotropic turbulence are available. In 
the following we will present results for these Reynolds 
numbers. All quantities are normalized by Kolmogorov 
scales, which means that they are non-dimensionalized by 
using the variables v and e. For example, all lengthscales 
are normalized by r/ = {v^/e)^^'^. This normalization 
allows to collapse the dissipation range of the different 
quantities if this range becomes independent of the vis- 
cosity. This is the case in the K41 phenomenology. In the 
presence of intermittency this is not the case anymore. 
It will however been shown in the following that also in 
that case the dissipation ranges of the different quantities 
nearly collapse in the present range of Reynolds numbers. 
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FIG. 1: Results for the energy spectrum computed by the 
EDQNM model. In the center plot we show the compensated 
spectrum. In the bottom figure the nonlinear transfer is plot- 
ted. All quantities are normalized by Kolmogorov scales. 



In Figure[l] the energy spectrum is shown for three dis- 
tinct Taylor-scale Reynolds numbers. We observe a clear 
^-5/3 power-law in the log-log representation. However, 
when showing the compensated spectra in log-lin repre- 
sentation it is observed that only at the highest Reynolds 
number a clear plateau can be discerned. At small k this 
plateau drops to zero, and at large k a viscous bottleneck 
is observed. 
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FIG. 2: Results for the energy spectrum computed from the 
multi-fractal description. In the center plot we show the com- 
pensated spectrum. In the bottom figure the nonlinear trans- 
fer is plotted. 



Since we are interested in second and third order quan- 
tities in the present work, we also show the nonlinear 
transfer. Again we observe that the asymptotic case, 
here indicated by a plateau around zero in between the 
negative and the positive lobe of the transfer spectrum, 
is only observed at the highest Reynolds number. 

The energy and transfer spectra computed from the 
multifractal description are shown in Fig. [51 The results 
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FIG. 3: The second order londitudinal structure function 
computed by EDQNM and the multifractal model. Straight 
black lines indicate powerlaw behavior proportional to and 



r2/3 



pensated according to K41 scaling. 



respectively. In the bottom figure the functions are com- 



are shown for relatively low Reynolds numbers (upto 
R\ = 400), since the numerical integration for higher 
values yielded extremely noisy results in the dissipation 
range. A bumpy large-scale behavior is observed in the 
compensated energy spectra, corresponding to the ad-hoc 
modeling of the large scales, as explained in appendix 
ICl This modeling also causes the relatively narrow neg- 
ative peak in the transfer spectrum. In the following, 
when presenting the structure functions, we do not need 
a smooth behavior for the large-scales and we will there- 
fore not use the ad-hoc modification of the large scales. 
In the dissipative and near-dissipative ranges, the spec- 
tral quantities (such as power spectrum and nonlinear 
transfer) obtained from EDQNM closures and the mul- 
tifractal formalism are very similar. In the dissipation 
range a viscous bottleneck is observed in both descrip- 
tions. 
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FIG. 4: The third order londitudinal structure function com- 
puted by EDQNM and the multifractal modeL Straight black 
hues indicate powerlaw behavior proportional to and 
respectively. In the bottom figure the functions are compen- 
sated according to K41 scaling. The asymptotic result 4/5 is 
indicated by a dashed line. 



Results for structure functions 



We used equations p6l) and ((20)) to compute the struc- 
ture functions from the energy spectra and transfer spec- 
tra shown in the previous section. The results for the 
second order structure function are shown in figure [31 We 
show the multifractal prediction in the same graph. In 
the log-log representation we clearly observe the smooth 

small scale behavior and the plateau proportional to 
the kinetic energy at large scales. In between a power-law 
dependence close to r^/"^ is observed. 

The multifractal prediction closely ressembles the 
EDQNM result in the dissipation range. The differences 
between the two models are more clearly visible in the 
compensated plot, where we observe that for R\ = 2500, 
the power-law dependence is clearly steeper than r^/"^. 
The largest difference is observed at large r. Indeed, the 
multifractal description does not take into account the 
shape of the velocity correlation at large r. This corre- 
lation should in a realistic flow smoothly tend to zero, 
but this effect is not taken into account in the formal- 



ism. Note that we prefer to show here the results with- 
out the ad-hoc modification proposed in the last section. 
We further observe that the structure function computed 
from EDQNM, as for the energy spectra, does not dis- 
play a clear plateau in the compensated representation 
for Rx < 25000. 

The results for the third-order structure functions are 
shown in figure |31 Again we clearly observe the smooth 
small-scale behavior proportional to . In this range the 
multifractal model closely follows the EDQNM results. 
For larger r a close to linear dependence and at large 
scales a decrease towards zero. Also here the multifractal 
formalism does not take into account the large scales. 
In the inertial range, at very large R\, the third order 
structure function should scale as 

DLLL{r) = -\er. (49) 
5 

it is observed that this is only reached at the highest 
R\ for EDQNM and only for a short range of scales. 
The multi-fractal results collapse with p9)) already at 
Rx = 380. 



Results for the velocity increment and 
derivative skewness 



In figure [5] top, we show the velocity increment skew- 
ness for different Reynolds numbers. In the K41 phe- 
nomenology, this quantity should give a constant value 
in the inertial range. It is observed that the fact that 
^LLhiy) tends to zero smoothly for large r results in a 
gently decreasing function, rather than a constant value. 

In the dissipation range all curves nearly collapse. 
Only the multifractal approach gives a slightly higher 
value than the rest, since the velocity derivative skew- 
ness is a function of Rx-, as will be shown later, in figure 
El At large scales, the multifractal result closely follows 
the high- Reynolds EDQNM result up to the cut-off of the 
multi-fractal result. 

In the center and top graph of figure [5l we compare 
the results also with experimental results. At Rx — 380, 
we compare with hot-wire measurements in an air-jet ex- 
27[. At small scales the experimental value is 



periment 

significantly larger than the theoretical results. At these 
scales the accuracy of the hot-wire probe decreases how- 
ever. In the inertial range the multi-fractal approach is 
very close to the experimental value. The EDQNM curve 
drops much faster to zero. Inhomogeneity and anisotropy 
of the experimental turbulent field could be behind this 
discrepancy. 

At Rx = 2500 a comparison is made with the veloc- 
ity increment skewness measured in the return-channel 
of the ONERA wind-tunnel in Modane. The Reynolds 
number obtained there is one of the highest measured 
in wind-tunnel turbulence. Unfortunately at large scales 
the third-order statistics are not fully converged so that 
no smooth curve is available there. However the gen- 
eral trend of the curves is quite similar at all scales. A 
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FIG. 5: Comparison of the longitudinal velocity increment 
skewness between EDQNM and the multifractal approach. 
In the center figure results are compared to air-jet results. In 
the bottom plot the results are compared to high-Reynolds 
number wind-tunnel experiments. The straight black line in- 
dicates a powerlaw proportional to r""'"* 



surprising fact is here the power-law that is observed in 
the inertial range of both the multifractal result and the 
EDQNM result. Indeed, in the multifractal approach this 
power-law is a signature of inertial-range intermittency 
and the model is developed to take this into account. 
In the EDQNM approach, however, this power-law is a 
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FIG. 6: The velocity derivative skewness as a function of the 
Reynolds numbers for EDQNM and MF. Also shown are DNS 
results from references [2^. [29|. 

transient effect, due to the finite-Reynolds number. In 
EDQNM this power-law vanishes thus at high Reynolds 
number. 

In figure [6] we show the Reynolds number dependency 
of the longitudinal velocity derivative skewness, as com- 
puted by equation We see that this quantity sat- 
urates for R\ > 100 at a value around 0.4. In the mul- 
tifractal approach this quantity follows a power-law of 
the Reynolds number (see Eqs. (|C7|) and (jC9p ) with 
an exponent around 0.13. Also shown arc the results of 
Direct Numerical Simulations [1^, [2^ . Note that in [2^ 
more numerical and experimental results are available in- 
cluding the experimental compilation by j30| . We chose 
however those which give the general trend. The DNS 
results show a slightly increasing trend from 0.4 to 0.6 
for a Reynolds number going from 10 to 1000. 



D. Influence of large-scale forcing 

To conclude this results section we address the influ- 
ence of a large scale forcing on the scaling of the velocity 
increment skewness. Indeed, in experiments of nearly 
isotropic turbulence we often consider a turbulence gen- 
erated by a grid, advected by a mean velocity. This cor- 
responds in the frame moving with the mean flow, to 
freely decaying turbulence. Direct numerical simulations 
of isotropic turbulence are often forced at the large scales 
in order to obtain a as high as possible Reynolds num- 
ber. The difference between the two types of turbulent 
flows is important. For example in 311] it was shown 
that the normalized dissipation rate is nearly twice as 
high in decaying turbulence as it is in forced turbulence. 
Also for the appearance of scaling ranges this difference 
can be important. The difference between decaying and 
forced turbulence in approaching the asymptotic form of 
the third-order structure function was reported in [s^ . 
In figure [7] we show how the inertial range scaling of the 
energy spectrum improves when considering statistically 
stationary forced turbulence at the same Reynolds num- 
ber. We observe the appearance of a large peak in the 
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FIG. 7: The influence of a large scale forcing on scaling are il- 
lustrated for the energy spectrum and the velocity increment 
skewness. The black-dotted lines correspond to forced tur- 
bulence at the same Reynolds number as the decaying cases 
considered. 



spectrum, corresponding to the forcing. We also observe 
that the compensated energy spectra display a clear scal- 
ing range, already at a Reynolds number of Rx = 380. 
The velocity-increment skewness of these forced calcu- 
lations shows however no clear plateau at this Reynolds 
number, but its inertial range behavior follows closely the 
R\ = 2500 decaying turbulence result. 



the multifractal model almost collapsed with the scaling- 
correction induced by the finiteness of the Reynolds num- 
ber in the EDQNM simulations. In particular at a 
Reynolds of R\ = 2500 the two corrections almost co- 
incide. This shows that at Reynolds numbers currently 
achievable in controlled experiments and simulations, in- 
termittency corrections to the skewness can not be dis- 
tinguished from low Reynolds number effects. An inter- 
esting perspective is to investigate to what extent inter- 
mittency corrections to higher-order quantities such as 
the flatness can be distinguished from Reynolds number 
effects (see e.g. Ref. [l^). This task is within the frame- 
work of closure-theory far from trivial and will be left for 
future work. 



Appendix A: Relation between second order 
structure function and the energy spectrum 

We will give here a detailed derivation of the relation 
between second order structure functions and the kinetic 
energy spectrum. Parts of this derivation can be found 
in different textbooks, but we think it is useful for the 
interested reader to give all the details in this work. 

Starting from p^ . homogeneity allows to write 



(Al) 
(A2) 



Using the inverse Fourier transform, and © we can relate 
this to the energy spectrum 



DllIt) = 2- 



UiUj 



,(fc)Mj(-fc)e*'=''dfc()A3) 



UiUj 



^M^E{k)e'''-^dk] (A4) 



ink 



Defining (f> the angle between k and r, we find 



r 



"^B^jik) (l-cos^, 



(A5) 



V. CONCLUSION 

In the present work we computed second and third or- 
der structure functions from EDQNM results. We com- 
pared these structure functions with results from the mul- 
tifractal formalism. It was shown that in the near dis- 
sipation range the different approaches give very similar 
results. It was shown that the appearance of clear scaling 
ranges is very slow for the structure functions as was also 
observed in previous work [l^, . 

The results for the velocity increment skewness were 
also compared to experimental results. It was shown that 
the intermittency correction to this quantity given by 



Also, in isotropic turbulence, the Reynolds stress tensor 
takes the form. 



UiUj 



^u^S^j, u^='^J E{k)dk, 



(A6) 



so that, introducing conveniently oriented spherical co- 
ordinates, we write, 

DLL{r) = 2^-2 [ '^^ ^(fc)e'^-"2^fc'sin(/)d,^dfc 



cos^ fl!))e''" '' sin 



dk. 

(A7) 
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The integral over (j) can be performed analytically by in- 
troducing C = cos (p and x — kr: 



(1 -cos^0) sm(t> e^'^'-^^'^c 



(1 - C^) e"'^dC 

L 

dx'^ J ix 



2(1 

4 



9^ \ sin a; 



dx'^ J X 
sin a; cos x 



yielding 



4 1 i?(fc) 



1 sin(fcr) — (fcr) cos(fcr) 
3 



(fcr)3 



(A8) 



Appendix B: Relation between third order structure 
function and the energy transfer spectrum 



The Fourier transform of ulu'j^ with respect to r is 

(Bl) 



UL{k)u\{-k) 



(B2) 



UL{k) J J UL{p)uLiq)S{-k - p - q)dpdq (B3) 
J J UL{k)uL{p)uL{q)5{k + p + q)dpdq (B4) 
■ jj Ui {k)uj {p)um {q)S{k +p + q)dpdq 



— 3 -^ijm\"'J 

Analogously we find 
So that 



3 ijjnK / 



(B5) 
(B6) 



(B7) 



DLLLir) = 1 (r„,„(fc) - 7^;„(fc)) e^'^-dfc. 

(B8) 

It is clear from the definitions p2|) and (|13p that 
(Tijm(fc) — T!*„(fc) is a third order tensor, symmetric 
in the indices j, m and solenoidal in the index i, so that 
its general form is 



Tijm{k) T^A^{k) — T{k)Pijrn{k) 



(B9) 



after multiplication of both sides by Pij,n{k) one finds 



r(fc) 



4fc2 

We substitute this in (IBS 



(BIO) 



D,,,ir) = 3^ / g^P.,„(fe)e--dfe. (BU) 
Defining the angle between k and r, we find that 



mjpHp..^ (fc) ^ 2k cos 0(1 - cos^ (j)) (B12) 



and thus 



DLLLir) = 6 /cos0(l -cos2(/))lA^e*'-dfc. (B13) 

Introducing again conveniently oriented spherical coordi- 
nates, we write this as 



DLLLir) = 6 / / cos 0(1 
Jo Jo 

27rfc^ sin (/)d0dfc 

T(fc) 



cos2 0)^^e^'=- X 



(B14) 



-3i 



sin 0d0dfc. 



cos (7!)(1 — cos (h)e x 



(B15) 



As for DLLir), the integral over can be performed an- 
alytically by introducing C. — cos0 and x — kr: 



cos(p (1 - cos^ (p) smcj) e*'"'"' 



C(l - C') e^^dC 



-2i 



2 sin X 6 cos x 6 sin x 



(B16) 



yielding 

DLLLir) 



, . , 3 (sin kr — kr cos kr) — [kr)'^ sin fcr 

lAT I 1. ( rC ) 7"; Qrv 



(fcr)5 



(B17) 



Appendix C: The multifractal description 

The multifractal formalism can be seen as a proba- 
bilistic interpretation of the averaged behavior of ve- 
locity structure functions. More precisely, for a scale 
r in the inertial range, using both the standard argu- 
ments of the multifractal formalism \9\ and the proba- 
bilistic formulation of Castaing 34], the velocity incre- 
ment 5uLix,r) — ULix + r) — ULix) can be represented 
by the product of two independent random variables, 
5uLix,r) = Pr^, with ^ a zero average Gaussian ran- 
dom variable of variance = {[SuLix, Lq)]'^), where 
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Lq is the integral length scale, and a stochastic vari- 
ance /3r — (xir) where the exponent h fluctuates itself 

according to the law Vrih) oc ijj^j . This gives 

a complete one-point probabilistic description (including 
structure functions and probability density functions) of 
the velocity increments in the inertial range given by an 
empirical function 'D[h). This function is both scale and 
Reynolds number independent, and is called the singular- 
ity spectrum in the inviscid limit. Paladin and Vulpiani 
[2JI and Nelkin [2^ then proposed a natural extension 
to the dissipative scales and the respective description 
of the velocity gradients. This adds to the description a 
Reynolds dependence through the fluctuating nature of 
the dissipative scale r]{h) = Lo{ne/Tl*)-^/'^''+^\ where 
T^e = cr^o / ^ Ei-nd TZ* = 52 a universal constant related to 
the Kolmogorov constant [13] . The relation between this 
Reynolds number and the Taylor-scale Reynolds number 



IS 



n* 



-R 



A- 



(CI) 



Let us remark that in a K41 framework the variable h = 
1/3 is unique and does not fluctuate and one recovers 
the classical Kolmogorov prediction rjj^ = r]{h = 1/3) = 

Lo(7^e/7^*)-3/4. 

Meneveau (ssj proposed an elegant interpolation for- 
mula between the inertial range and the far dissipative 
range. Following these works reference [13] proposed a 
probabilistic formulation of velocity increments that cov- 
ers the entire range of scales. The expressions for and 
Vr{h) in this description are 



and 



Vr{h) 



Tl(h) 



-|(l-/i)/2 



(C2) 



Z(r) 



1 



r,{h) 



{V{h)-l)/2 ' 



(C3) 



where Z{r) 

rh 



is a normalization factor such that 

Given the parameters of the flow, namely Lq, the large 
scale variance and TZe, this description requires one 
additional free parameter TZ* and a parameter function 
2?(/i) that can be measured from empirical data. We 
will take TZ* = 52 and a parabolic approximation for the 



singularity spectrum 'D{h) 



1 



2C2 



The so-called 



intermittency coefficient has been estimated from data 
to be C2 = 0.025 [13] • The remaining parameter ci is 
chosen such that, in the inertial range, the third order 
structure function (|5ui]'^) is proportional to the scale r. 
This gives ci = 



§C2. 



The proposed description has been shown to accurately 
describe the symmetric part of the velocity increments 
probability density functions and even order structure 
functions. In particular, using Eqs. (jC2p and (|C3[) . even 
order structure functions are given by the following inte- 
gral 



/3^'^Pr{h)dh , (C4) 



with (^^'') — o'^'^-^sil. Furthermore, one can show [1? 
that the mean dissipation 



lim 



LqTZ* ^ ' 



is independent on the Reynolds number. Since the Gaus- 
sian noise ^ is independent on the fluctuating exponent 
/i, odd order structure functions vanish. Some modifica- 
tions of the noise ^ have been proposed in Refs. [13] and 
[13] in order to take into account a non-zero skewness. 
Nevertheless, without any additional free parameters, the 
Karman-Howarth-Kolmogorov equation 



-(e)r + Qv 



d{{5uLf) 
dr 



(C6) 



gives a prediction for the third order moment {(dui^)^), 
and thus the skewness, knowing only the second order 
one. This study was carried out and compared against 
experimental data in Ref. [27| . As a final remark, we 
would like to add that such an approach allows to give a 
prediction for the third order moment of velocity gradi- 
ents, i.e. 



2(0) 



[2/1 - 1 - V{h)] ( -V 1 dh + F 



L 



(C7) 



where is a negligible additive term, coming from the 
Taylor's development of the normalization factor Z{r), 
and given in [27|- Using Eq. (jC7p . a steepest-descent 
calculation shows that the skewness of the derivatives 
behaves as a power law of the Reynolds number, i.e. 
-S'(O) - 7^^-l with 



X = mm 

h 

3 . 

mm 

2 h 



2{h - 2) + \ - V(h) 



h+l 
2{h- 1) + 1 



Vih) 



h + l 



(C8) 



Using a quadratic approximation for the parameter func- 
tion V and C2 — 0.025, one gets 



' S{0) - TZ 



0.067 



(C9) 
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These result were already obtained by Nelkin [2^ 
using a different, although related, approach, based 
on the asymptotically exact relationship {{dxu)^) — 
-M{dluf)- 

Finally, we need to discuss the behaviour of the mul- 
tifractal description at scales of the order of the inte- 
gral length scale L^. At this stage, the proposed for- 
malism (Eqs. (IC2[) and (jCSp ) is valid only in the limit 
r ^ Lq, such that the integrals correspond to their 
steepest-descent values. The multifractal description de- 
scribes hereby an asymptotic state in which the influ- 
ence of the decay of the turbulence or the energy-input 
mechanism is not taken into account. The inertial range 
starts directly at the integral scale. Second order struc- 
ture functions do therefore not smoothly tend to a con- 
stant value around Lq and the third-order structure func- 
tion does not tend to zero. The present work is not 
devoted to an extension of the description to take into 
account explicitly the large scales. However, in order to 
compute the energy and transfer spectra from structure 
functions through the relations (|18|) and ((22|) we need a 
smooth behavior around the large scales for the integrals 



to converge. For this reason we propose an ad-hoc mod- 
ification of the multifractal description. Firstly we use 
Pr = 1 (i.e., independent on h and non fluctuating) and 
Vih) uniform for r > L. This is equivalent to a Gaus- 
sian modeling of the velocity increments, with a scale 
independent variance. Unfortunately, this description is 
not continuous (because of the dissipative corrections). 
Moreover, one has to extend the validity of the Karman- 
Howarth-Kolmogorov equation (|C6I) itself in order to get 
a realistic behaviour of the third order moment at large 
scales. In order to get a continous and differentiable in 
scale description of the velocity increments, we propose 
to replace the scale r entering Eqs. (|C2[) and (|C3|) by 
the ersatz r' = Lq tanh(r/Lo)- The hyperbolic tangent 
allows a smooth transition to r' = r for r ^ Lq and 
r' = Lq for r ^ Lq. Still, this does not fix the unrealistic 
behavior of Eq. (jC6[) at large scales. Thus, the predic- 
tion of {[6u{x,r')]'^) obtained from Eq. (|C6P using the 
ersatz r' is furthermore multiplied by a large-scale cutoff 
of the form exp(— r'^/(2LQ)). We stress again that this 
approach is completely ad-hoc and only used in order to 
allow computation of wave-number spectra. 



[20] 



A. Kolmogorov, Dokl. Akad. Nauk. SSSR 30, 301 (1941). 
L. Richardson, Proc. Roy. See. London, Series A 110, 
709-737 (1926). 

J. Lumley, Phys. Fluids 4, 206 (1992). 

L. Landau and E. Lifschitz, Fluid Mechanics, Addison- 

Wesley, 1959. 

A. Kolmogorov, J. Fluid Mech. 13, 82 (1962). 

A. Oboukhov, J. Fluid Mech. 13, 77 (1962). 

C. V. Atta and W. Chen, J. Fluid Mech. 44, 145 (1970). 

F. Anselmet, Y. Gagne, E. J. Hopfinger, and R. A. An- 
tonia, J. Fluid Mech. 140, 63 (1984). 

U. Frisch, Turbulence, the legacy of A.N. Kolmogorov, 

Cambridge University Press, 1995. 

R. Kraichnan, J. Fluid Mech. 5, 497-543 (1959). 

R. Kraichnan, Phys. Fluids 8, 575 (1965). 

R. Kraichnan, J. Fluid Mech. 47, 513 (1971). 

Y. Kaneda, J. Fluid. Mech. 107, 131 - 145 (1981). 

S. Orszag, J. Fluid Mech. 41, 363 (1970). 

W. Bos and J.-P. Bertoglio, Phys. Fluids 18, 031706 

(2006). 

W. Bos, H. Touil, and J.-P. Bertoglio, Phys. Fluids IT, 
125108 (2005). 

J. Qian, Phys. Rev. E 55, 337 (1997). 
J. Qian, Phys. Rev. E 60, 3409 (1999). 

G. K. Batchelor, The theory of homogeneous turbulence, 
Cambridge University Press, 1953. 

A. Monin and A. Yaglom, Statistical Fluid Mechanics II, 
MIT Press, 1975. 



[21 

[22; 

[23; 

[24; 

[25' 
[26; 

[27 

[28; 

[29 
[30 

[31 

[32; 

[33 

[34 

[35; 



J. Mathieu and J. Scott, An introduction to turbulent 

flow, Cambridge University Press, 2000. 

M. Lesieur, Turbulence tn fluids, Kluwer Dordrecht, 

1990. 

P. Sagaut and C. Cambon, Homogeneous Turbulence Dy- 
namics, Cambridge University Press, 2008. 

G. Paladin and A. Vulpiani, Phys. Rev. A 35, 1971 
(1987). 

M. Nelkin, Phys. Rev. A 42, 7226 (1990). 

H. Kahalerras, Y. Malecot, Y. Gagne, and B. Castaing, 
Phys. Fluids 10, 910 (1998). 

L. Chevillard, B. Castaing, E. Leveque, and A. Arneodo, 
Phys. D 218, 77 (2006). 

T. Ishihara, Y. Kaneda, M. Yokokawa, K. Itakura, and 

A. Uno., J. Fluid Mech. 592, 335 (2007). 
R. M. Kerr, J. Fluid Mech. 153, 31 (1985). 

K. R. Sreenivasan and R. A. Antonia, Annu. Rev. Fluid 
Mech. 29, 435-472 (1997). 

W. Bos, L. Shao, and J.-P. Bertoglio, Phys. Fluids 19, 
045101 (2007). 

Y. Kaneda, J. Yoshino, and T. Ishihara, J. Phys. Soc. 
Japan 77, 064401 (2008). 

L. Chevillard, B. Castaing, and E. Leveque, Eur. Phys. 
J. B 45, 561 (2005). 

B. Castaing, Y. Gagne, and E. Hopfinger, Phys. D 46, 
177 (1990). 

C. Meneveau, Phys. Rev. E 54, 3657 (1996). 



